Models of environment and T\ relaxation in Josephson Charge Qubits 



Lara Faoro, 1 Joakim Bcrgli, 2 '0 Boris L. Altshuler, 2, 3 and Y. M. Galperin 4, 5,6 

1 Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway 08854, New Jersey, USA 
2 Physics Department, Princeton University, Princeton, New Jersey 0854-4, USA 
3 NEC Laboratories America, Inc., 4 Independence Way, Princeton, NJ 08540, USA 
^Department of Physics, University of Oslo, PO Box 1048 Blindern, 0316 Oslo, Norway 
5 A. F. Ioffe Physico- Technical Institute of Russian Academy of Sciences, 194021 St. Petersburg, Russia 
"Argonne National Laboratory, 9700 S. Cass av., IL 60439, USA 
(Dated: February 2, 2008) 

A theoretical interpretation of the recent experiments of Astafiev et. al. on the Ti-relaxation rate 
in Josephson Charge Qubits is proposed. The experimentally observed reproducible nonmonotonic 
dependence of Ti on the splitting Ej of the qubit levels suggests further specification of the previously 
proposed models of the background charge noise. From our point of view the most promising is the 
"Andreev fiuctuator" model of the noise. In this model the fiuctuator is a Cooper pair that tunnels 
from a superconductor and occupies a pair of localized electronic states. Within this model one can 
naturally explain both the average linear T\{Ej) dependence and the irregular fluctuations. 

PACS numbers: 



Proposals to implement qubits using superconducting 
nanocircuits have undergone an amazing development 
during the last years 0, Q, S S 0| ■ In Josephson Charge 
Qubit (JCQ) information is encoded in the charge states 
of a Cooper pair box. The JCQ is manipulated by tuning 
gate voltage and magnetic flux. Both time resolved co- 
herent oscillations in single and coupled JCQ have been 
recently observed 0, El ■ Although decoherence is a se- 
vere limitation to the performances of these devices the 
dominant source of noise is yet to be identified. 

A significant step towards a characterization of the en- 
vironment in a JCQ has been recently made by Astafiev 
et.al. The experimental set-up consists of a Cooper 
pair box connected to a reservoir through a tunnel junc- 
tion of SQUID geometry with Josephson energy Ej 
pierced by an external magnetic field. Provided that 
E c ^> Ej ^> T (where E c , Ej and T are correspond- 
ingly charging energy, Josephson energy and tempera- 
ture, ks = ft = 1)- only two charge states |0) and |1) arc 
relevant and the Hamiltonian of the box reads: 



H q = 



Ej 



(1) 



E c (l — CgVg/e), Cg is the gate capacitance, 



where SE, 

In the rotated basis {|+), |— )} the Hamiltonian reads 



g is the gate voltage and e denotes the electron charge. 



E 



Pz 



(2) 



where E = y 8E 2 C + Ej and 6 = arctan(£j /SE C ). One 
can distinguish the off degeneracy working points (9 ~ 
and 6E C S> Ej) and the degeneracy one (6 = tt/2 and 
5E C = 0). Astafiev et.al. |7J measured the energy re- 
laxation rate T\ of the JCQ in a wide range of parame- 
ters. Two main features have been observed: (*) Linear 



increase of Y\ with Ej at large Ej, and (**) Small non- 
monotonous fluctuations in the ri(£',/)-function on this 
linear background. 

We do not believe that the existing experimental infor- 
mation is sufficient to identify a unique interpretation. 
However, it substantially reduces the range of possibil- 
ities. In this Letter we show that some models which 
have been used to study dephasing in JCQ can not ex- 
plain these features. We propose a model where all of 
them appear naturally. 

Many different mechanisms can be responsible for de- 
coherence in JCQ. We will consider three models, all 
based on the idea that the oxide layer close to some 
metallic reservoir, like one of the leads or gates or Cooper 
pair box itself, is disordered and thus hosts trapping cen- 
ters, i.e. localized states for the electrons. Model I: The 




Model I Model II Model III 

FIG. 1: The three models 

reservoir is a normal metal, and electrons can tunnel from 
any state below the chemical potential to an unoccupied 
trap above the chemical potential or from an occupied 
trap to an extended state above the chemical potential. 
Model II: There is no tunneling between the reservoir 
and the traps, but an electron in an occupied trap be- 
low the chemical potential can be excited into an empty 
trap. Model III: The reservoir is a superconductor and 
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a Cooper pair can split in such a way that each of the 
electrons end up at an empty trap (Andreev fluctuator). 
At the end of the paper we discuss why we do not be- 
lieve that dephasing through phonon or photon modes 
can describe the experimental results. 

To describe the qubit interacting with the environment 
one has to supplement Eq. by the Hamiltonians of the 
environment, He, the interaction Hi, and the tunneling 
Ht- Regardless to its particular features an environment 
is coupled to a JCQ through the charge degree of free- 
dom: a z = p z cos 9 — p x sin 9. Since only the term p x in 
the interaction Hamiltonian will change the state of the 
qubit, the relaxation rate Ti should be proportional to 
sin 2 9 if the qubit-environment coupling is weak. 

The charge of the qubit affects the environment in two 
ways. First, the electrostatic interaction shifts the energy 
of the localized states. Second, the amplitude t of the 
tunneling between the trap and the reservoir (another 
trap) depends on the qubits charge: t(o~ z ) = to + ta z . 

The total Hamiltonian can be written as 

H = H q + Hi + H t + H E . (3) 
For the environment Hamiltonian we write 

a k 

n E ~ Zj ° 

a 

where c a , (c£) destroys (creates) an electron in trap a 
and fk, (fl) destroys (creates) an electron in the reser- 
voir. (We assume that the superconductor in Model III 
is always in the ground state). 

Ht describes the tunneling in the absence of the qubit 

H^=t J2i<ifk+c a ft), 

k j a 

h[ h) = h E <&<*, 

The coupling of the qubit with the environment is gov- 
erned by the Hamiltonian 

Hj= (vJ^clca+^H^J a z (4) 

Let us now discuss processes where the qubit, initially 
prepared in the excited state, releases the energy by ex- 
citing the environment (Ti-proccsscs). When the cou- 
pling is weak the relaxation rate can be derived by using 
the Fermi Golden Rule. If originally the qubit is in the 
excited state and the bath occupies a state \i) with prob- 
ability pi the decay rate to the ground state, Tj, is 

T i = 2n*£p°\(+,f\Te- i ti H >W\-,i)\ 2 . (5) 

ij 



Here |/) is the final state of the bath. If the 
qubit is prepared in the ground state the transi- 
tion rate to the excited state is at thermal equilib- 
rium Tj = e~ E / T Ti. The total relaxation rate is then 
T 1 = T L + T T = (1 + e~ E/T )T i . In the limit E > T we 
have Ti w r x . 

When using Eq. JSJ the exponential has to be ex- 
panded to the necessary order. Notice that the second 
term in Eq. J3J that describes the change in the tun- 
neling amplitude will both flip the qubit (remember that 
o~z = Pz cos 9 — p x sin 9) and change the occupation of the 
trap. Thus, it contributes to Ti already in the first order: 

rf 5 = 2irt 2 g(E) sin 2 9. (6) 

Here g(u>) is the density of states of excitations. The elec- 
trostatic interaction term, '^o~ z c*c, does not change the 
occupation of the localized state, so it contributes only 
in the second order. Assuming that (togo) 2 /E « g(E), 
where go is the density of states in the metal, we can 
write this contribution as 

if 5 = ^v 2 [to + 2tcos9] 2 g(E) sin 2 9 (7) 
E* 

For Model II, the contribution J7J) does not appear as 
long as all traps are coupled equally to the qubit. Then 
moving one electron from one trap to another will not 
change the electrostatic potential. Accordingly for Model 
lithe v 2 should be interpreted as (v 2 ) averaged over some 
scatter of v. 

Let us calculate g(oj) for Models I- III. Consider the 
density of localized state: 

v{e) = v + 5v{e), (8) 

where v is the average value of i/(e) (we assume that v is e- 
independent). The random deviations Sv(e) are assumed 
to be small, 5v(e) v and only short-range correlated: 

{5v{t)5v{(')) = A8(e - e') (9) 

For the density of states (g{u))} averaged over different 
realizations of the random distribution of trap energies, 
we have 

g(w)W = go [ dev(e), (<?H) (/) = 2Dg u 

g{u) {II) = I" fev{e)v{e-w) t (g(u)}W=D 2 w 
Jo 

<?M (//7) = f 'de[v(e)v(u> - e) (10) 
Jo 

+ u(-e)u(e - u)], {g{u)) {III) = 2v 2 u. 

go denotes the density of states in the metal and we ne- 
glect its energy dependence. 
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In each of the three models (g(oj)) is a linear function 
of the frequency ui. 

This statement may become incorrect when Coulomb 
interaction between trapped electrons is taken into ac- 
count. For example in Model II the tunneling ampli- 
tude depends exponentially on the distance between the 
traps. Only traps which are close in space can exchange 
charge. For such pairs the Coulomb interaction between 
the traps modifies the density of states to (see Ref. 0) 



<sH> 



(a) 



u 2 (e 



(ii) 



lu). where e c is the energy of the 



Coulomb interaction between two electrons that occupy 



the two traps, e^ 1 ^ 1 can be estimated as e^ 11 ^ 



2 /r u 




where rt is the typical distance between the traps, i.e. 
typical tunneling length 

As to the Model III the Coulomb repulsion leads to 
(g{tu)) {III} = 2v 2 (lu - e { c III} ). When estimating e { c IH) 
we need to take into account the screening provided by 
the superconductor ||: each electron trapped near the 
superconductor creates an image charge and thus forms a 
dipole moment of the order of ert- The distance between 
the two traps is determined by the coherence length £ of 



the superconductor. Therefore e 



(in) 



It is informative to compare the Coulomb energies 
^ with the superconducting gap A ~ Hvf/S,, where 

VF 



FIG. 2: The relaxation rate as function of frequency for Model 
III. The curve shows one typical realization of the disorder 
while the solid line is the average over realizations. The 
dashed line shows (Ti(lu) 2 ). 



straight lines represent the result of an ensemble averag- 
ing. If the fluctuations are rapid one can instead average 
over suitable frequency windows for a single sample. 

The second order term can also give rise to fluctuations 
in the relaxation rate. Let us focus on Model I at the op- 
timal point (cos 9 = 0). The main source of fluctuations 
is the w 2 -term and for this it follows immediately from 



I I ^ u i - / 3/ iB Liic v -teini anu iui tins it lunuws 

is the Fermi velocity Since e 2 /hv F « 1, we Eqs q and j^j that the correlator is 



have e { c n '/A ~ £/r t , and e ( c ln '/A ~ (r t /£) 2 . It is nat 
ural to assume that r t <C £. Therefore A <C e^ 7 ' and 
When considering JQC one is interested 



A > e 



(in) 



(12) 



Ej <C A. We conclude that lu <C e 



(ii) 



and 

(Hi) 



thus (g(u>)) ( - 11 ' 1 = const. At the same time lu e c 
{gito))^ 111 ^ is determined by Eq (|TU)l and the "Andreev 
fluctuators" can lead to a linear dependence of Ti on Ej. 

Returning now to the equations JHJl and (0 we see that 
the first order contribution is directly proportional to the 
density of states. The second order term, because of the 
E in the denominator, does not give a linearly increasing 
relaxation rate even if the density of states is linear. 

We conclude that to get a linear rate from a linear 
density of states we need some term in the Hamiltonian 
that gives a contribution to first order. 

So far we only considered the average relaxation rate, 
let us now turn to the fluctuations. 

Consider the two-point correlator: 



Thus, also in Model I there will be fluctuations, but the 
peaks will have a different shape, and the correlations are 
long range. Also, in Models II, III the amplitude of the 
oscillations increases with increasing E, whereas in Model 
I it decreases since it only has contributions to second 
order. This could be a way to distinguish the different 
models. For comparison with Figure|2is shown in Figure 



( 5 ( w ) 5 (^'))i /7/) 



2Agl min(wV), 

2Av 2 min(w, w') + A 2 lo8(lo - o/jjjl) 
%Av 2 min(w,a/) + AA 2 loS(uj - lu'). 




Note that from Eqs. (fTUfl and (fTTfl it follows that in Model 
I g(Lu) is a monotonous function of w, whereas Models II, 
III dead to non-monotonous fluctuations with short range 
correlations. 

Figure shows for Model III at the optimal point the 
relaxation rate as function of frequency for a particular 
realization of the position in energy of the traps. The 
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FIG. 3: The relaxation rate as function of frequency for Model 
I. The curve shows one typical realization of the disorder while 
the solid line is the average over realizations. The dashed line 
shows (ri(aj) 2 ). 



13 Ti for Model Hot the case where v 2 t^jl 2 
units of E. 
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Above wc took the correlation of the levels to be a 
true 5-function. In reality, the levels will be broadened 
by relaxational processes that go beyond our models. To 
take this into account we can use instead for the corre- 
lator some (5-like function (5 CT (e) with characteristic width 
<7. For example, one can think about Gaussian, 5^ = 
e -t' 72er 2 j^p2^ a or Lorentzian, 5 { a L) = a/n{e 2 + a 2 ), cor- 
relations. The (5-functions in Eq (jl 1|> are then replaced 
by the function 5& which again is a (5-like function; 
in the Gaussian case <5"i = <5 ( ' < 2 and for Lorentzian, 

v2cr 



J 2cr 
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Note that both phonon and photon radiation could 
cause a linear frequency dependence of the relaxation rate 
in two dimensions due to their linear dispersion. How- 
ever we do not believe that they can be responsible for 
the observed resonances. In this case a peak in the Ti 
as function of E follows directly from a resonance in the 
density of states g{E). Let us estimate frequency of such 
a resonance assuming that the resonant structure in the 
density of states arise from quantization ofphonon levels 
in a confined geometry. According to Ref. [jj one such res- 
onance was at a frequency of 30 GHz. The sound velocity 
is 10 m/s and corresponds to a wavelength of 30 nm. 
While not impossibly small, this appears to be smaller 
than the typical sizes of > 100 nm of the structures in 
the samples used. On the other hand the possibility that 
there could be a coupling to a standing photon mode in 
the experimental cavity looks more likely. A similar ar- 
gument but using the speed of light gives us a wavelength 
of 1 cm, which is of the right order of magnitude. Only 
two samples where measured Q with slightly different 
resonant frequencies (20 and 30 GHz), but this could be 
caused by different position or size of the sample. In 
view of the fact that the cavity contains the sample and 
mount as absorbing material and that no special care was 
taken to create a high Q cavity it appears unlikely that 
such a sharp resonance line would be created. This could 
be tested by introducing some absorbing material into 
the cavity and see if the resonant peaks will disappear. 
An alternative way to discriminate between a phonon or 
photon resonance peak and one created by a resonant 
fluctuator would be to thermocycle a given sample. If 
the resonance is caused by some fluctuator, the latter 
probably would be rearranged by the heating, and thus 
the peak positions would shift. 

In summary, we have argued that dephasing by 
phonons or photons is unlikely to explain the experimen- 
tal results although they can not be ruled out conclu- 
sively. A more likely explanation is some resonant fluc- 
tuator model. We have discussed three such models, and 
all of them depend on the effect of the state of the qubit 
affecting the barrier height to reproduce the linear de- 
pendence of the relaxation rate on E. 

We think that Model III (Andreev fluctuators) is the 
most promising for the following reasons. Models II, III 



allow for rapid, nonmonotonous oscillations of the Ti for 
large E, whereas Model I to first order only will show 
steplike monotonous increase of Fx. To second order 
there are nonmonotonous oscillations also for Model / but 
they have a different shape. Model II is less likely than 
Model III because the Coulomb interaction most likely 
changes the density of states to constant for the relevant 
range of energies. 

To experimentally determine the coupling constants 
we suggest the following: If one probes the same en- 
ergy E at different working points (by changing both 
6E C and Ej) there should to first order be collapse of 
the data points if one plots Fi/sin 2 # as a function of 

(2) 

E, whereas the terms with cosd in T\ will cause some 
deviation. In particular, it seems instructive to plot 
[T(E,6)/T(E,tt/2) - l]/cos0 = 4(v/E) 2 (t /t + cos6) as 
function of cos#. From this one could extract the ratios 
v/E and t /t. 

We proposed thcrmocycling as an experimental check 
for the presence of fluctuators and introduction of some 
absorbing material in the cavity to rule out photon reso- 
nances. 
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